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Using a first-principle pseudo-potential approach, we have investigated the Schottky barrier heights 
of abrupt Al/Ge, Al/GaAs, Al/AlAs, and Al/ZnSe (fOO) junctions, and their dependence on the 
semiconductor chemical composition and surface termination. A model based on linear-response 
theory is developed, which provides a simple, yet accurate description of the barrier-height variations 
with the chemical composition of the semiconductor. The larger barrier values found for the anion- 
than for the cation-terminated surfaces are explained in terms of the screened charge of the polar 
semiconductor surface and its image charge at the metal surface. Atomic scale computations show 
how the classical image charge concept, valid for charges placed at large distances from the metal, 
extends to distances shorter than the decay length of the metal-induced-gap states. 

PACS numbers: 73.30.+y, 73.40.Ns 



I. INTRODUCTION 

Metal/semiconductor (MS) interfaces have been the 
focus of extensive theoretical and experimental studies 
for several decadesiS*^ To date, however, we are still 
far from a complete understanding of the factors which 
control the Schottky barrier height (SBH) at these in- 
terfaces. In recent years, new research activities have 
been developed in the area of band engineering at MS 
interfaces 5 * 6 * 7 *^^*^ and on the properties of metal/wide- 
gap-semiconductor contactsiii*i2ii£ These developments 
have stimulated renewed interest in some basic issues 
concerning Schottky barriers, and in particular in the 
mechanisms that control the SBH dependence on bulk- 
semiconductor and interface-specific characteristics. 

The problem of Schottky barrier formation has been 
traditionally addressed by studying the dependence of 
the SBH on the metal used in the junction^ Early stud- 
ies suggested a Schottky-Mott behavior controlled by the 
metal work function for highly ionic or wide-gap semicon- 
ductors, and a weak dependence on the metal type and 
on the junction fabrication method for the most cova- 
lent semiconductors, such as Si or GaAs^i^ The latter 
trend was generally attributed to various Fermi level pin- 
ning mechanisms, such as pinning by metal-induced-gap 
states (MIGS)— at an intrinsic charge neutrality level of 
the semiconducto r 16 : 17 or pinning by native defect states 
of the semiconductor at some extrinsic gap levelii&iS Fur- 
thermore, a correlation between Schottky barriers and 
heterojunction band offsets was observed experimentally 
for a number of systems^ and similarly ascribed to 
Fermi level pinning at a bulk reference level. Finally, 
the effect of the semiconductor ionicity on the SBH trend 
with the metal work function was examined in pioneering 



self-consistent studies of jellium/semiconductor contacts, 
and the trend could also be generally understood in terms 
of MIGS properties of the semiconductor. 21 

More recent experiments on metal contacts to covalent 
semiconductors, however, have revealed a much weaker 
electronic pinning than was previously believed; 2 In par- 
ticular, there have been reports of considerable changes 
in metal/Si and metal/GaAs SBH's obtained by alter- 
ing the structural properties and/or the chemical com- 
position of the interface A 2 ^ The conclusion that the 
SBH does depend most generally on the microscopic 
atomic structure of the interface has been reached by 
many authors, both on experimentalSii^SiLSiiiii 2 ^ and 
theoretical 11 : 23 ' 24 : 25 ' 26 : 27 : 28 grounds. While opening a 
promising line of research on Schottky barrier engineer- 
ing, these observations complicate seriously the search 
for simple models of Schottky barrier formation, since 
the inclusion of the interface atomic structure seems un- 
avoidable. 

Given the complexity and variety of the atomic struc- 
ture at metal/semiconductor contacts, it seems unlikely 
that a simple unified model could emerge and entirely 
cover the various facets of Schottky barrier formation. 
Conversely, a systematic investigation of the problem 
starting from abrupt, defect free interfaces, and progres- 
sively introducing perturbations at the interface could 
help identifying relevant microscopic mechanisms and 
provide a firmer basis for modeling Schottky barrier prop- 
erties. Progress in computational physics have made 
possible accurate ab initio calculations of the electronic 
structure of MS contacts, and the complexity of the sys- 
tems which can be examined is steadily increasing; this 
type of computations can provide the means to carry out 
such an investigation and probe the correlation between 
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microscopic atomic structures and SBH's. The present 
study is a first step in this direction. 

In this article, we study from first principles the depen- 
dence of the SBH on selected bulk and surface character- 
istics of the semiconductor, for a given metal. Specifi- 
cally, we examine abrupt Al/X (100) junctions, where 
X = (Ge, GaAs, AlAs, and ZnSe) are lattice matched 
semiconductors of increasing ionicity, and we investi- 
gate the microscopic mechanisms responsible for the SBH 
changes with the semiconductor chemical composition 
and surface termination (cation or anion). A model based 
on a linear-response-theory scheme is then developed, 
which explains our ab initio results and SBH trends ob- 
served experimentally. 

II. METHOD OF CALCULATION 

We have carried out ab initio calculations, within the 
Local-Density Approximation (LDA) to Density Functio- 
nal Theory (DFT), using the pseudopotential method^ 
We used norm-conserving scalar-relativistic Troullier- 
Martins pseudopotentials2°*2I in the Kleinman-Bylander 
non-local formal and the exchange-correlation functional 
of Ceperley and Alder^ The electronic states were ex- 
panded on a plane-wave basis set using a kinetic energy 
cut-off of 20 Ry. We used supercells containing 7 Al layers 
and 13 semiconductor layers (7 + 13 supercell) to model 
defect-free Al/X (100) junctions. In section lTV Bl we em- 
ployed larger supercells (7 + 21) to investigate the screen- 
ing of substitutional charges placed in the junctions, and 
to compute the parameters (D s , 5 S ) necessary to model 
this screening. All supercell calculations were performed 
with a (2, 6, 6) Monkhorst-Pack fc-point gridi^i 

We considered ideally abrupt epitaxial junctions and 
neglected atomic relaxation at the interfaces. The ef- 
fect of atomic relaxation at the Al/GaAs (100) and 
Al/ZnSe (100) interfaces has been examined in Refs. 11 
and |35l Atomic relaxation decreases (increases) the p- 
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(a) (b) 

FIG. 1: (a) Epitaxial alignment of Al on the (100) surface 
of zinc-blende semiconductors verifying the lattice-matching 
condition osai = a/y2. (b) Atomic structure of the abrupt 
Al/semiconductor (100) interface. The semiconductor surface 
is terminated either by an anion (I) or a cation (II) plane. 



type Schottky barriers of the abrupt Al/GaAs (Al/ZnSe) 
junctions by 0-0.1 e\4^ (0.1-0.2 eVilsifc), and has no in- 
fluence on the SBH ordering of the anion- and cation- 
terminated interfaces. 

The epitaxial alignment of Al on the (100) surfaces 
of the four semiconductors under study is illustrated in 
Fig. nja). This type of alignment corresponds to the 
lattice- matching condition: a ax — a /y2> where a is the 
semiconductor lattice parameter. The Al [100] direction 
is parallel to the semiconductor [100] axis, and the whole 
Al fee lattice is rotated by 45° about its [100] axis with 
respect to the semiconductor substrate. Experimentally 
- and also in our calculations — Ge, GaAs, AlAs, and 
ZnSe are lattice-matched semiconductors, and the equi- 
librium lattice constant of Al is slightly larger (1%) than 
a/y/2. This results in a small compressive strain in the 
Al in-plane lattice parameter, which is accommodated 
by an elongation (~ 3%) of the Al overlayer, assuming 
pseudomorphic conditions. For the semiconductor lat- 
tice parameter, we used the theoretical value a = 5.55 A 
(a ox P- = 5.65 A). The metal-semiconductor interlayer dis- 
tance d at the junction was taken as the average between 
the (100) interlayer spacings in the semiconductor and in 
the (strained) Al bulk parts, i.e., d = 1.72 A. The po- 
lar Al/X (100) junction offers two inequivalent abrupt 
interfaces, either with anion- or cation-terminated semi- 
conductor surface, which we both considered in our study 
[see Fig.^b)]. In what follows, we will refer to the anion- 
and to the cation-terminated interface as interface I and 
II, respectively. 

To evaluate the p-type SBH, <p p , we used the same 
approach as in previous studies >2Mai 

p = AV + AE p , (1) 

where AV is the electrostatic-potential lineup at the in- 
terface and AE p is the difference between the Fermi level 
in the metal and the valence-band maximum (VBM) in 
the semiconductor, each measured with respect to the 
average electrostatic potential in the corresponding crys- 
tal. The band-structure term AE p is characteristic of 
the individual bulk crystals forming the junction. This 
term was computed using the Kohn-Sham (KS) eigen- 
values of standard bulk band-structure calculations. The 
potential lineup AV contains all interface-specific contri- 
butions to <j) p and was obtained — via Poisson's equation 
and using a macroscopic average technique^S^ — from 
the self-consistent supercell charge density. 

For a meaningful comparison of our calculated SBH's 
with experiment, AE p should include quasiparticle and 
spin-orbit corrections. The spin-orbit correction is simply 
+-% 2 ) where A so is the total spin-orbit splitting at the 
semiconductor valence-band maximum, which was taken 
from experiment. For a metal, in principle, the exact KS 
Fermi energy and the quasiparticle Fermi energy must 
coincide at zero temperature. 37 Furthermore, LDA cal- 
culations for the work functions of various Al surfaces 
performed with the same method and the same pseu- 
dopotentials as in the present study — and neglecting 
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many-body corrections on the Al Fermi energy — yielded 
values which agree with the experimental data to within 
a few tenths of meVi£& In the present study, therefore, for 
the metal Fermi energy we just used the LDA result. The 
corrected band term is thus AE p — AE^ S — AE qp — , 
where AEp^ 3 is the KS band term and Ai? qp is the dif- 
ference between the quasiparticle and KS semiconductor 
VBM energies. 

For the quasiparticle corrections, we used the results of 
GW calculations taken from the literature. For Al/GaAs 
we used the correction AE^ aAs = —0.36 eV evaluated 
by Charlesworth et a/., 25 who employed for the refer- 
ence LDA calculations the same exchange-correlation 
potential as we do. For Ge (AlAs) we used the cor- 
rection for GaAs, and the difference between the Ge 
(AlAs) and GaAs corrections evaluated in Ref. l39l i.e., 

A£7° e - AS^aAs = + q Q 9 e y (A£AlAs _ A £GaAs = 

—0.11 eV)4& The quasiparticle corrections to the band 
structure of ZnSe have been evaluated in Ref. l4ll As 
the LDA bandgap in our calculations and in Ref. |4l] are 
different, due to the different pseudopotentials e mplo yed, 
we took the valence-band-edge correction of Ref. |4l] and 
scaled it by the ratio of the difference between the LDA 
and GW bandgap in the two calculations. The resulting 
estimate for AE qp is —0.50 eV. Using the experimental 
spin-orbit splittings A A = 0.30, 0.34, 0.28, and 0.43 eV 
for X — Ge, GaAs, AlAs, and ZnSe^ the total correc- 
tions are 0.17, 0.25, 0.36, and 0.36 eV, respectively^ 
The numerical uncertainty on the absolute value of the 
SBH's is estimated as ~ 0.1 eV for Al/Ge, Al/GaAs, and 
Al/AlAs, and as ~ 0.2 eV for Al/ZnSe. For a given inter- 
face geometry, however, the relative barrier values (A<fi p 
in Table HT|) are considerably more accurate, i.e, have an 
estimated numerical accuracy of ~ 50 meV. 

III. RESULTS FOR THE SCHOTTKY BARRIER 
HEIGHTS 

The calculated SBH's for the abrupt Al/X (100) inter- 
faces, including many-body and spin-orbit corrections, 
are given in Table [I] We observe a systematic differ- 
ence between the type-I and type-II interfaces: the p-type 
SBH is always higher for the type-I (anion-terminated) 
interface. This difference increases with increasing semi- 
conductor ionicity. Our theoretical results are compared 
with experimental SBH values in Fig. [5] For the Al/Gc, 
Al/GaAs, and Al/AlAs systems the experimental ranges 
correspond to data obtained by transport measurements. 
In the case of Al/GaAs, photoemission measurements — 
performed at low metal coverage — give rise to a wider 
range of SBH values^ but the scattering in the data de- 
creases significantly when thick metallic overlayers are 
deposited and the barriers are measured by transport 
techniques. For Al/ZnSe we are not aware of any trans- 
port data and we used photoemission results. 

In the case of Al/Ge, no SBH measurement has been 
performed, to our knowledge, on the (lOO)-oriented inter- 



TABLE I: Estimated quasiparticle and spin-orbit correc- 
tions to </>p DA for different semiconductors. The calculated 
Al/X (100) SBH's including these corrections are shown in 
the last two columns. All numbers are in eV. 



Semiconductor 
X 


Estimated 
correction 


<t> P 

I 


II 


Ge 


0.17 


0.21 




GaAs 


0.25 


0.86 


0.76 


AlAs 


0.36 


1.45 


1.16 


ZnSe 


0.36 


2.18 


1.82 



face. In Fig. |21we have thus used the existing transport 
data 4 ^ on Al/n-Ge (111) junctions {<j> n = 0.52-0.61 eV), 
together with the Ge experimental bandgap at room tem- 
perature, E^ e = 0.66 eV4& The resulting barrier heights 
4> p = 0.05-0.14 eV compare reasonably well with our cal- 
culated value of 0.21 eV. In the case of Al/GaAs (100) 
the transport measurements give values of 4> p between 
0.58 eV and 0.76 e V 8 i 9 i 10 i 47 i 48 i 49 i 50 i 51 i 52 i 53 i 5 V 5 i 56 (we 
used ijGaAs _ ^ 42 eV to estimate the p-type barrier 
heights from measurements performed on Al/n-GaAs 
junctions). This is in relatively good agreement with 
our calculated SBH of 0.76 eV for the Ga-terminated in- 
terface, and still consistent with our value of 0.86 eV 
for the As-terminated interface. Concerning the effect of 
the GaAs-surface stoichiometry on the measured SBH, 
we note that different conclusions have been reached by 
different groups. Some studies, including Refs. and 
I49L have found a small (~ 0.1 eV) difference between 
the SBH's measured in junctions fabricated on As-rich 
and on Ga-rich surfaces (As-rich leading to higher 4> p , 
consistent with our results), while other studies, such as 
Refs. l5fl and I53L found no difference. 

For the Al/AlAs system, Ref. Hj] reports SBH values 
for Al/n-AlAs (100) ranging from 0.85 to 0.94 eV for var- 
ious reconstructions of the semiconductor surface, while 
somewhat higher values, 0.95 and 1.01 eV, have been 
given in Refs. |5^ and l56l respectively. Using the exper- 
imental bandgap £" A1As = 2.16 eV, the resulting range 
is 4> p = 1.15-1.31 eV, in good agreement with the calcu- 
lated SBH for the Al-terminated AlAs surface (1.16 eV), 
and somewhat smaller than the value we find for the 
As-terminated surface (1.45 eV). The Al/n-ZnSe (100) 
SBH has been investigated in Refs. 0, 0, and |53 for 
different reconstructions of the ZnSe (100) surface. Very 
similar values have been reported for the c(2 x 2) and 
2x1 reconstructions, namely <j) p = 2.12-2.15 eV n i 12 i 57 
and <j> p = 2.11-2.15 eVjii*^ respectively, while a lower 
SBH, <pp = 1.91 eV, has been measured for the lxl 
reconstructioniii These values are in between our values 
of 1.82 eV for the Zn-terminated and of 2.18 eV for the 
Sc-tcrminated interface. 

The general agreement between theory and experiment 
in Fig. |21 indicates that our calculations for ideal MS 
structures (abrupt interfaces with no atomic relaxation) 
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FIG. 2: Schottky barrier heights (SBH's) at Al/X (100) con- 
tacts, X — Ge, GaAs, AlAs, and ZnSe. The circles show the 
calculated SBH's for the ideal anion- and cation-terminated 
junctions, corrected for quasiparticle and spin-orbit effects. 
The shaded regions show the ranges of experimental values 
(see text). 



capture the general trend of the SBH with the chemi- 
cal composition of the semiconductor. We note that for 
Al/ZnSe, the inclusion of the appropriate reconstruction 
and relaxation brings the theoretical results in very close 
agreement with the experimental values^ 



IV. INTERPRETATION AND MODELS 

A. General trend with the semiconductor chemical 
composition 

Experimentally a correlation was found between Schot- 
tky barriers and heterojunction band offsets^ A large 
number of MS contacts and semiconductor heterojunc- 
tions were shown to verify, within 0.4 eV, the transitivity 
relationship: 



<t> p [M/Si] - <f> p [M/S 2 ] = AE YBO [S 1 /S 2 



(2) 



where M is a metal (such as Al or Au; in general nei- 
ther a highly reactive nor a transition metal 5 ^), S\ and 
S% are two semiconductors, and AEVbo [S1/S2] is the 
corresponding valence-band offset (VBO). This correla- 
tion was most often observed for MS junctions used in 
transport measurements, i.e., which had been annealed 
for fabrication of the contacts. The experimental data in 
Fig. |21 are in general agreement with the above empirical 
transitivity rule. 

We note that the transitivity rule, as formulated in 
Eq. (J2J , disregards any dependence of the SBH on the mi- 
croscopic interface structure, and cannot therefore give 
a complete account of the theoretical results in Fig. [2] 
Also, recent theoretical and experimental studies have 
shown that the band offset at heterovalent semiconduc- 
tor heterojunctions depends critically on the orientation 



and other microscopic details of the interface ^SMi The 
right-hand side of Eq. J2J is thus ill defined, in general, 
for heterovalent semiconductors. 

In this section we concentrate on the average SBH 
i'p = 3 \$p + ^p 1 ) °f t ne abrupt, defect-free type-I and 
II interfaces, and propose a model for its variation with 
the semiconductor chemical composition, derived from 
an atomic-scale approach. We show that this variation 
is controlled essentially by the same bulk mechanism 
that determines band offsets at non-polar, defect-free 
semiconductor heterojunctions^ 9 - The splitting A(j> p = 
(fip — (fip 1 due to the semiconductor-surface termination 
will be the focus of the next section. 

Similarly to the SBH, the VBO may be written as: 
ASvbo = AE V + Ay, where AE V is the difference be- 
tween the VBM energies of the two semiconductors, each 
measured relative to the mean electrostatic potential in 
the corresponding crystal, and AV is the electrostatic po- 
tential lineup at the interface. Since the band-structure 
terms AE V and AE p [in Eq. JTJ] are differences between 
bulk values of the individual crystals forming the junc- 
tion, they verify by definition the transitivity relationship 
in Eq. All non transitive contributions are contained 
thus in the potential lineup terms AV. 

In the case of semiconductor heterojunctions a linear- 
response-theory (LRT) approach, which focuses on AV 
and treats the interface as a perturbation with respect to 
a bulk reference system, has provided an accurate general 
description of band-offset trends. 59 6 - Based on this ap- 
proach and comparision with fully self-consistent ab ini- 
tio calculations, it has been shown, in particular, that in 
the case of defect-free, isovalent lattice-matched semicon- 
ductor heterojunctions, AV is determined by the prop- 
erties of the bulk constituents (as opposed to interface- 
specific features, such as interface orientation or interface 
abruptness). Specifically, if S% and 5*2 are the two semi- 
conductors, with anion (cation) species di (ci) in Si and 
&1 (C2) in 52, the potential lineup is given within LRT 
by»S 



AV [5i/5 2 ] = 



2ire 2 
3fl 



r 2 [An a (r) + An c (r)}dr, (3) 



where the integration is over the whole space, Q is the vol- 
ume of the bulk unit cell, and An a (An c ) is the electronic 
charge density induced by a single anion (cation) substi- 
tution ai — ► a 2 (ci — > c 2 ) in the bulk semiconductor 5iiS 
Based on this LRT approach, it has also been shown that 
in the case of heterovalent lattice matched semiconduc- 
tors, Eq. ® also applies in the specific case of defect-free 
interfaces with the non-polar (110) orientation^ 

Using a similar linear-response scheme for AV, we 
show in Appendix f£] that the average SBH <j> p can be 
described by the following model: 



jmod 
» p 



4> p {Al/(X) (100)] + A£ V bo[(A)/A (110)]. (4) 



The first term on the right-hand side of Eq. (0J is the 
SBH at the (100) interface between Al and the group-IV 
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FIG. 3: Comparison of the average SBH 0p DA at the 
Al/X (100) I and II interfaces with the model predictions, 
Eq. Q) (left panel) [Eq. (right panel)]. The horizontal 
bar shows the average SBH and the small dot indicates the 
model result, i.e., the sum of the SBH at the A\/{X) (100) 
[Al/Ge (100)] junction (gray circles [dashed line]) and the 
VBO at the (X)/X (110) [Ge/X (110)] interface (double ar- 
rows). The quasiparticle and spin-orbit corrections are not 
included; these contributions trivially verify Eqs. (0J and 0. 



virtual crystal, denoted (X), which is obtained by aver- 
aging the anion and cation pseudopotentials of the III-V 
or II- VI compound X (X = GaAs, AlAs, ZnSe). The 
second term is the VBO of the non-polar (X)/X (110) 
heterojunction. 

The basic approximation to derive Eq. Q is to con- 
struct the charge densities of the Al/X I and II junctions 
(and hence their average lineup) starting from the refer- 
ence Al/ (X) system, by adding a linear superposition of 
the charge densities induced in the virtual crystal (X) by 
single anion and cation substitutions that transform (X) 
into X. The Al/ {X) (100) junction is an optimal refer- 
ence system in this context, which minimizes the devia- 
tions of 4> p from 4> I p° d in Eq. (@}; these deviations vanish 
to the first order in the ionic substitutions which trans- 
form Al/(X) into the Al/X I and II junctions. 

It is also possible to use as a reference system another 
Al/group-IV (100) junction, whose density is sufficiently 
close to the average density of the Al/X I and II junc- 
tions. For instance, one may use Al/Ge as a common 
reference system and obtain (Appendix |A"|> : 



/mod' 
P 



b p [Al/Ge (100)] 



The deviations of 



A£ V Bo[Ge/V (110)]. (5) 

ftp from 4> p nod in Eq. © include, in 
this first order correction in the substitutions. 

The latter correction can be identified with the dipole 
induced in the reference Al/Ge (100) junction by isova- 
lent Ge — » (X) substitutions performed within the first 
one to three Ge atomic layers closest to the interface (see 
Appendix^); such a dipolar term is generally small for 
isovalent substitutions (~ 0.1 eV or less, see Ref . l2Sj) . and 
will be neglected here. 

In Fig. we compare graphically the model predic- 
tions, Eq. (0J and Eq. J5J, with the calculated average 



SBH of the Al/X I and II interfaces. The (110) VBO's 
have been computed using supercells containing 8 planes 
of each semiconductor in the ideal (unrelaxed) lattice- 
matched geometry. The same energy cutoffs and fc-points 
grids have been used as in the calculations of the Schot- 
tky barriers. The SBH of the Al/{X) (100) junctions 
have been obtained using the same parameters as for the 
Al/X (100) I and II interfaces. The results in Fig. 
show that Eq. (@J and Eq. © provide a fairly accurate 
(±0.15 eV) description of the average SBH (p p . We note 
that the SBH's at the Al/(X) junctions are all small, due 
to the small bandgaps of the virtual crystals (< 0.4 eV), 
and similar to the LDA SBH at the Al/Ge (100) interface 
(0.04 eV). The results in Fig.EJshow that the (110) VBO 
— a bulk-related quantity in our calculations — controls 
the general increase of the barriers from the group- IV to 
the III-V and to the II- VI semiconductors. 



B. Effect of surface termination 



We will show here that the difference A.<p p = <fi p — (ft^ 
due to the semiconductor-surface termination in Fig. [3 
and in particular the fact that the SBH is systematically 
higher for the anion than for the cation termination — 
can be understood in terms of surface-charge and image- 
charge effects. The mechanism is illustrated in Fig. 0] 

With respect to the Al/ (X) interface, the ionic charge 
distributions of the interfaces I and II are obtained by 
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FIG. 4: (a) Planar average of the difference between the ionic 
charge densities of the anion- (cation-)terminated Al/X (100) 
interface and the XI/ (X) (100) interface; a = 1 for the semi- 
conductors X — GaAs, AlAs, and a = 2 for X = ZnSe. (b) 
Macroscopic average of the ionic charge density difference, (c) 
Positive (negative) potential difference established at the in- 
terface I (II) by a positive (negative) surface charge and its 
image charge at the metal surface. 
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substituting an anion (charge +er) on each anionic site 
and a cation (charge —a) on each cationic site, as indi- 
cated in Fig. Ufa). For the III-V and II- VI compounds 
we have a = 1 and a — 2, respectively, in units of 
charge per unit-cell surface in the (100) plane. The ar- 
rows in Fig. Ufa) represent opposite delta functions on 
each anionic and cationic (100) plane, corresponding to 
the planar average of the ionic point-charge density. The 
macroscopic average 5 ' 9 '- of this ionic charge distribution 
is represented in Fig. E{b). In the bulk semiconductor, 
the macroscopic average eliminates the atomic-scale os- 
cillations of the planar charge density; at the interface, 
however, a positive (negative) charge density subsists in 
the junction I (II). This macroscopic charge has a den- 
sity p = ±2a/a and extends over a distance a/4 be- 
tween the last Al plane and the first semiconductor plane. 
It is therefore equivalent to a surface charge of density 
a = ±a/2. 

Within a classical macroscopic description, a plane of 
charge in a semiconductor is screened by the dielectric 
constant of the host material. Furthermore, in the 
presence of a metal the screened surface charge is neu- 
tralized by an image charge induced at the metal surface, 
and a potential difference is thus established between the 
two charges [see Fig.^Jc)]. If a is the density of surface 
charge in the semiconductor, x the position of the plane 
of surface charge, and Xi the position of the metal surface 
or image charge, the potential difference obtained from 
classical electrostatics is: 



AU(x, a) = Aire 2 — (x - Xi). 



(6) 



As can be seen from Fig. 0Jc), in the junction I such a 
dipole lowers the average potential energy in the semi- 
conductor with respect to its value in the metal, increas- 
ing the SBH <p p ; conversely, in the junction II the dipole 
raises the average potential in the semiconductor, de- 
creasing <p p . The mechanism illustrated in Fig. 0] thus 
provides a qualitative explanation for the difference be- 
tween the SBH's of the interfaces I and II. Of course, 
the classical limit given by Eq. © is expected to be cor- 
rect only for a charge placed at a large distance from the 
metal. As we will see below, however, closer to the metal 
the above type of description may still be used provided 
the inhomogeneous nature of the screening near the metal 
is taken into account. 

To check that the mechanism in Fig. 0] can indeed ac- 
count for the SBH differences A0 p , we have calculated 
the changes in the lineup (and hence in the SBH) in- 
duced by surface charges of varying magnitude, placed 
on the semiconductor plane closest to the metal in the 
Al/X I and II junctions. At the interface I (II), a surface 
charge of density — |<r| (+|er|) was introduced by replac- 
ing the anion A (cation C) of the semiconductor layer 
adjacent to the metal surface by a virtual ion (Ai-z.Cz.) 
((Ci-z A-)). The resulting changes AU in the SBH ob- 
tained from the ab initio calculations are shown in Fig. [5] 
The negative (positive) surface charge decreases (in- 



creases) the p-type SBH of the anion-terminated (cation- 
terminated) interfaces, consistent with the screened sur- 
face charge and image charge description in Fig.Efc). We 
also note that, consistent with the latter description, the 
bare monopole is replaced by an interface dipole in the 
multipole expansion of the total (electronic plus ionic) 
charge disturbance. 

The macroscopic average of the difference between the 
ionic potentials in the junctions II and I, in Fig. 01 a), is 
equivalent to a surface charge er = — 1 at the interface 
for the III-V semiconductors and a — — 2 for the II- VI 
semiconductors. Therefore, focusing on the effect of the 
macroscopic charges only and to the first order in the 
perturbation, the modification of the SBH in the junc- 
tions I for a = —1 (—2) should be equal to the difference 
(ftp 1 — (ftp for the III-V (II- VI) semiconductors. Similarly, 
the change of the SBH in the junctions II induced by a 
surface charge a = +1 (+2) should be equal to (ft p — (ft p l . 
Our ab initio results in Fig. [S] show, however, that the 
responses AU of the two interfaces are not linear when 
\o~\ > 0.5 and differ in magnitude. Therefore, we take the 
average AU between the potential differences induced in 
the junctions I and II as our estimate for the difference 
A(j) p . The results are shown in Table ITT1 For Al/ZnSe 
we reported the calculated SBH changes for a = ±2 (not 
shown in Fig.[SJ). The average values AU are seen to de- 
scribe well the calculated difference A0 p and also the in- 
crease of A(j} p when the semiconductor ionicity increases. 

Although this supports the surface charge and image 
charge picture in Fig.QJc), Eq. needs to be revisited 
for charges placed close to the metal surface. For ex- 
ample, in the case of a test charge a = 0.1 on the first 
semiconductor layer, we obtain from Eq. © a potential 
difference AU = 82 meV, using for the distance x — Xi the 
value d/2, where d = 1.72 A is the interplanar distance 
at the interface (see Fig. ^| , and for eoo the theoretical 
dielectric constant of GaAs, e^ aAs = 12.4^ This result 
is more than 5 times larger than the ab initio result for 



Anion-terminated 




Cation-terminated 



1.0 



■ Al/ZnSe 
a Al/AlAs 
• Al/GaAs 



-0.5 0.0 0.5 

a (electrons/S) 

FIG. 5: Schottky barrier modification AU induced by a bare 
surface charge a per unit-cell surface S = a 2 /2 on the first 
semiconductor layer closest to the metal in the Al/X (100) I 
and II junctions (see insets; the same symbols for the atomic 
layers are used as in Fig. 0. 
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TABLE II: Comparison of the average SBH change AU in- 
duced by surface charges a — ±1 (GaAs, AlAs) and a — ±2 
(ZnSe) at the interface, with the difference Acj> p = (f> p — 
(ftp 1 between the SBH of the anion- and cation-terminated 
Al/X (100) junctions (see text). The last column shows the 
results of the model Eq. I1H . All energies are in eV. 



X 


A0 P 


|AU| 


AU 


D s 


Ss 


A(f> p aod 






I II 




(eV" 1 !" 2 ) 


(A) 


Eq. (nj 


GaAs 


0.10 


0.18 0.06 


0.12 


0.051 


2.5 


0.10 


AlAs 


0.29 


0.27 0.19 


0.23 


0.060 


2.0 


0.13 


ZnSe 


0.36 


0.68 0.24 


0.46 


0.041 


1.8 


0.40 



the potential difference obtained for Al/GaAs and shown 
in Fig. 

In Ref. |2(] we observed a somewhat similar behavior in 
the case of local dipoles inserted in the Al/GaAs (100) 
junction. In the latter case, ab initio calculations were 
performed to determine the change in the lineup Au in- 
duced by a dipole layer (i.e., two test charges +<r and 
—a placed on two adjacent cation-anion planes) intro- 
duced at various distances x from the metal, within the 
semiconductor. In this numerical experiment, the bare 
dipole perturbation is Aw& = Aite 2 o~l, where I is the sep- 
aration between the charged planes, and from the com- 
puted Am we could directly measure the effective dipole 
screening e^P = Aub / Au as a function of the dipole po- 
sition x in the junction. This screening was found to be 
strongly inhomogeneous and to increase exponentially as 
the dipole was approaching the metal surface. This was 
attributed to the MIGS tails and their high polarizability 
in the interface region. We also proposed a model for e^P 
which proved very accurate to describe the SBH changes 
Au(x, a) in the linear-response regime (i.e., to the first 
order in ct)>2& 



ejj? (x) w eoo + A-Ke 2 D s (E F , x) 5 S 



(7) 



Here D s (Ep, x) is the MIGS surface density of states at 
the Fermi energy and at the position x of the dipole, and 
S s is the decay length of the MIGS. The model Eq. is 
also consistent with earlier MIGS-based model descrip- 
tions of Schottky barrier properties iSLSiS 

In order to predict, in general, the effect of the surface 
termination in MS junctions, we would like to develop a 
model for AU that takes into account the inhomogeneous 
nature of the electronic screening in the MIGS region 
and that is consistent with our previous results on the 
effect of the dipole layers on the SBH. In particular, this 
model should be consistent with the fact that, in the 
linear-response regime (small |cr|), the sum of the SBH 
modifications induced separately by two charges +a and 
—a separated by a small distance I, AU(x — 1/2, a) and 
— AU(x + l/2, a), respectively, must be equal to the SBH 



modification induced by the corresponding dipole: 

AU(x - 1/2, a) - AU(x + 1/2, a) = (8) 

e effW 

Expanding the left-hand side of Eq. (JHJ to the first order 
in I, we obtain the differential equation: 



dAU 

dx 



Ane 2 cr 
e eff» 



(9) 



With our expression for e^(x) in Eq. J7J and a 
surface density of states that decays exponentially, 64 
D s (Ep, x) = D s (Ep, 0) exp(— x/S s ), the solution of 
Eq. © with the boundary condition AU(xq) = is: 



AU(x, a) = Aire 2 



x -5 s log 



4s ( x o) 



dip 

£ cff 



(x) 



, (10) 



where x is the position of the surface charge. We note 
that for large values of x, we recover the classical limit 
given by Eq ©, with Xi = x + 5 S log[ej! 1 f p (xo)/e 00 ]. 

With the exception of xo, all parameters necessary to 
evaluate AU(x, a) in Eq. i|10|) can be obtained straight- 
forwardly from ab initio calculations performed either 
for the bulk semiconductor (eoo) or for the unperturbed 
A\/X junction (the MIGS-related parameters). In Ta- 
ble [HJ we have reported our calculated values for the 
MIGS parameters D s = D s (Ep, 0) and 5 S . These quan- 
tities were obtained from the calculated macroscopic av- 
erage of the local density of states, N(E, x), as D s = 
J Q N(Ep, x)dx and 8 S — j Q xN(Ep, x)dx, where 
the origin (x = 0) was taken as the midpoint between 
the last Al and the first semiconductor plane, and oo 
indicates a position well inside the semiconductor (the 
center of the semiconductor slab in the supercell) where 
the MIGS vanish. As the values of D s and S s are slightly 
different for the interfaces I and II, we reported in Ta- 
ble [H] the average between the values calculated for the 
two interfaces^ 

In order to obtain an estimate for xq, and also to test 
the model in Eq. (|10fl . we have investigated ab initio the 
spatial dependence of AU in the linear-response regime 
by introducing a small test surface charge a — ±0.05 
in the As-terminated Al/GaAs junction at different dis- 
tances from the interface. This was done by replacing 
single layers of As (Ga) ions by virtual (Aso.95Sio.05) 
((Gao.95Sio.05)) anions (cations). As an example, we 
show in Fig.Elthe ab initio results for the charge density 
and potential induced by such a test charge on the sixth 
semiconductor layer from the metal. The macroscopic 
averages of the ionic, electronic, and total charge densi- 
ties are displayed in Fig.||Ja). We have used a Gaussian 
filter function with full width at half maximum a/2 for 
the macroscopic average. This allows one, in particular, 
to distinguish the image charge contribution to the total 
charge density, close to the Al surface. The macroscopic 



8 



n - 



- Image 


„ (a) - 


charge t 






- 
- 


\ V 

, x 10 " 


- Al As 


\J' Ga~ 




. ^ a 95Sio.05 


At/ 








- o = u.uo 


-5 


5 10 




x (A) 




FIG. 6: (a) Macroscopic average of the electronic (thin 
solid line) and ionic (dotted line) charge densities induced 
by a plane of (Gao.95Sio.05) virtual ions in the As-terminated 
Al/GaAs (100) junction. A Gaussian filter function was used 
for the macroscopic average. The thick solid line is the sum of 
the electronic and ionic densities, scaled by a factor of 10. (b) 
Macroscopic average of the corresponding induced total elec- 
trostatic potential. The resulting potential difference AU is 
also indicated, (c) Schottky barrier modification \AU\ (filled 
squares) obtained for a surface charge \a\ — 0.05 as a function 
of its position within the semiconductor, in the As-terminated 
Al/GaAs junction. The symbols give the results of the self- 
consistent calculations. The solid line corresponds to the pre- 
diction of Eq. il l (H with xq — 0.6 A. The atomic positions are 
indicated using the same symbols for the atoms as in Fig. 
The calculations were done in a 7 + 21 supercell. 



average of the induced total electrostatic potential is dis- 
played in Fig. E{b). The corresponding potential differ- 
ence is AU = 0.24 eV. 

In Fig. HJc) we have plotted the discontinuity \AU\ 
induced by the test charge as a function of its position 
in the Al/GaAs junction. With the theoretical dielectric 
constant of GaAs (e§f As = 12.4) and the calculated val- 
ues of D s and S s given in Table ITT1 the best fit of Eq. I|1(J|) 
to the data in Fig.|SJc) is obtained with xq = 0.6 A. The 
model results obtained from Eq. I|1U|I using xq — 0.6 A 
have been reported in Fig. EJc), and compare well with 
the results of the self-consistent calculations as a function 
of the distance. 

Having a reasonable estimate for xq we may now use 
Eq. (|10|l to obtain also an estimate for the difference be- 
tween the SBH's of the interfaces I and II. As we have 
seen before, this difference may be evaluated as the po- 
tential change induced, to the first order, by a surface 
charge a = 1 (2) per unit-cell surface on the first plane 
of the III-V (II- VI) semiconductor, i.e., at a position 
x = d/2 with d = 1.72 A (see Fig. [TJ. The resulting 
estimate A0™ od for the difference A(j> p is thus 



mod 



47T6 



d/2 -x -S s log 



e cfr P <>o) 



dip 
£ cff 



(d/2) 



(11) 



The results obtained with this model are displayed in 
Table HU To evaluate A<^ od , we have used the theo- 



retical value of the semiconductor dielectric constant e^, 
(X = GaAs, AlAs, and ZnSe)^& the surface density of 
states D s and MIGS decay length <5 S given in Table |H] 
for the Al/GaAs, Al/AlAs and Al/ZnSe junctions, and 
xq = 0.6 A for all systems. The model results, in Ta- 
ble [nj yield the correct trend and order of magnitude for 
the difference between the SBH's. We note that decreas- 
ing Xq increases A</>™ od , but does not affect the trend. 
The same conclusion applies when a possible change in 
xq, from GaAs to AlAs and to ZnSe, is taken into account 
by scaling the xq value obtained for GaAs with the ratio 
of the MIGS decay length, i.e., 5* /Sf aAs & 

Although the present model approach provides a con- 
sistent picture of the effect of selected perturbations at 
MS interfaces, we would like to caution the reader that 
our model description applies to unrelaxed interfaces. In 
this connection it should be pointed out that a LRT 
model (based on dynamical effective charges) is avail- 
able in the literature for quantitative predictions of the 
atomic relaxation contribution to the SBH. 27 



V. CONCLUSION 

Using a first-principle pseudopotential approach, we 
have investigated the Schottky barrier heights of 
abrupt lattice-matched Al/Ge, Al/GaAs, Al/AlAs, and 
Al/ZnSe (100) junctions, and their dependence on the 
semiconductor chemical composition and surface termi- 
nation. The sensitivity of the SBH to microscopic inter- 
face features reveals the limits of the currently accepted 
semiempirical models of Schottky barriers. Such model 
theories generally neglect the effects of the microscopic 
interfacial morphology. This is due in part to the com- 
plexity of the actual atomic structure of most MS con- 
tacts, and also to the relatively limited information avail- 
able on the atomic-scale geometry of buried interfaces. 

Based on our ab initio studies, we have derived mod- 
els which explicitly include the effects of the interface 
atomic structure in the case of some prototype defect- 
free, lattice-matched geometries. These models retain, 
within specific ranges of applicability, the same accuracy 
as the ab initio calculations. They show, in particular, 
that while the variation of the average SBH of the abrupt, 
anion- and cation-terminated Al/semiconductor (100) in- 
terfaces can be explained mainly in terms of the bulk 
properties of the semiconductors, the difference between 
the barrier heights of the anion- and cation-terminated 
interfaces results from a microscopic dipole generated by 
the screened charge of the polar semiconductor surface 
and its image charge at the metal surface. Our atomic- 
scale computations have also allowed us to show how the 
classical image charge concept, valid for charges placed 
at large distances from the metal, can be extended to 
distances shorter than the decay length of the metal- 
induced-gap states. 
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APPENDIX A: MODEL FOR THE AVERAGE 
SBH OF THE ANION- AND 
CATION-TERMINATED INTERFACES 

In order to explain the behavior of the average SBH 
(j) p with the semiconductor composition, we have ex- 
tended to MS contacts a linear-response-theory approach 
commonly used to study band offsets at semiconductor 
heterojunctionsp£2i§S The present analysis is also an ex- 
tension to heterovalent semiconductors of an approach 
outlined in Ref. |23 to model the Schottky-barrier changes 
with the alloy composition in Al/Gai_ x Al x As junctions. 

We are interested in the average potential lineup AV = 
\ (AV 1 + AV 11 ) , where AV I(II) is the potential lineup 
at the interface I (II), for which we want to establish the 
following model: 

AV m AV [A1/(Y) (100)] + AV [(Y) /X (110)] . (Al) 

The first term on the right-hand side of Eq. (|A1JI is the 
potential lineup at the reference (100) junction between 
Al and a group- IV (real or virtual) crystal (Y) (e.g., the 
virtual crystal (X) or Ge), having a charge density close 
to the average charge density of the Al/X I and II junc- 
tions. The second term is the lineup at the non-polar 
(110) interface between the group-IV crystal (Y) and the 
semiconductor X. To derive Eq. IjAlfl . we write the self- 
consistent electrostatic potential at the Al/X (100) I (II) 
junction as 

V l( - n \r) = V (r)+Vl {Il \r), (A2) 

where Vo(r) is the electrostatic potential at the Al/(Y) 
(100) junction. The average lineup AV can be expressed, 
according to Eq. (£2j, as AV = AVq+AVx, where AV = 
AV [Al/(Y) (100)] and AV X is the lineup of the potential 

= -[*?(»") + *?(»•)]• (A3) 

The potential V^ 1V> is the self-consistent electrostatic po- 
tential induced in the Al/ (Y) (100) junction by the ionic 
substitutions, (Y) — ► anion and (Y) — ► cation, performed 
in the semiconductor, which transform the A1/(Y) sys- 
tem into the type-I (type-II) Al/X system. Thus V/ and 
Vj 11 have long-range contributions associated with each 
heterovalent anion and cation substitution in the group- 
IV crystal. These long-range terms cancel out in the av- 
erage in Eq. I|A3() . since each anion (cation) substitution 



in Vl is compensated by a cation (anion) substitution 
associated with the same site in V/ 1 . The average poten- 
tial V\ has therefore a well defined macroscopic average 
in the semiconductor, which is equal to AV\, since V\{r) 
vanishes in the metal. 

One may thus evaluate AVi using a perturbative ap- 
proach neglecting inter-site interactions in the ionic sub- 
stitutions, because of the short-ranged nature of the po- 
tentials associated with each individual site. Within this 
approximation, V\ is given by i) the superposition of the 
potentials induced by isolated anion and cation substitu- 
tions in the bulk crystal (Y) plus it) a correction due to 
the deviations from the bulk response for substitutions 
performed near the MS interface. 

By construction, the potential lineup obtained from 
i) is orientation independent and equal to the potential 
lineup AV [{Y)/X (110)] at the non-polar (Y)/X (110) 
interface, built from a superposition of the same isolated 
charge-density responses on one side of the X (110) ho- 
mojunction. Furthermore, previous ab initio and LRT 
studies of semiconductor heterojunctions have shown 
that the deviation of the lineup AV (110) from the tran- 
sitivity rule, and hence from AV (110), is typically less 
than 0.1 eV in IV-IV/III-V junctions, and of the order 
of 0.1 eV in IV-IV/II-VI j unctions -^ffi We may there- 
fore replace AV [{Y)/X (110)] by AV [(Y)/X (110)] to 
obtain the contribution from i) to AVi. 

The correction to the lineup induced by ii) is given, to 
the first order in the substitutions, by AV corr . — J^jdj, 
where di = Aire 2 J dx xArii(x) is the dipole, and 7ij(x) 
the charge density, induced by the (Y) — > anion or (Y) — > 
cation layer substitution within the i th atomic plane from 
the interface in the A1/(Y) junction. In practice, the 
dipoles di vanish beyond the 3 rd to 4 th atomic plane 
from the junction, and AT4 orr . is generally of the or- 
der of 0.1 eV for isovalent substitutions^ Furthermore, 
when the Al/ (X) junction is used as a reference sys- 
tem, AV^orr. exactly vanishes, because the corrections 
are opposite in the I and II junctions and cancel out 
in the average leading to A V\ . The correction AV COYI . is 
therefore bound to be small (^0.1 eV) when the refer- 
ence system is an A1/(Y) (100) junction with a density 
close to the average density of the Al/X I and II junc- 
tions. We will thus neglect this correction, which leads 
to AVi ~ AV [{Y)/X (110)], and hence to Eq. pHjl. 

In Table II I II the average potential lineup at the 
Al/X (100) I and II interfaces is compared to the predic- 
tions of the model, Eq. i|Al|) . obtained with (V) = (X) 

and with (Y) = Ge. The agreement between AV m ° 

[AV m ° d ] and the calculated AV is quite good, the dis- 
crepancy being 2% [8%] or less when the Al/(X) (100) 
[Al/Ge (100)] junction is used as a reference system. 
Introducing the band energies in Eq. JM}, we obtain 
Eq. Q [Eq. ©] with (Y) = {X) [(Y) = Ge]. 
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TABLE III: Comparison of the average potential lineup AV at the Al/X (100) I and II interfaces with the predictions of the 

model, AV « AV[A1/(F) (100)] + AV[{Y)/X (110)] [Eq. glj], using (Y) = (X) (AF mod ) and (Y) = Ge (AF mod '). All 
numbers are in eV. 



X 


AV 


AV 




AF mod 


Ay 




AF mod ' 






Al/(X) (100) 


(X)/X (110) 




Al/Ge (100) 


Ge/X (110) 




GaAs 


-2.34 


-2.18 


-0.12 


-2.30 


-2.05 


-0.28 


-2.33 


AlAs 


-2.01 


-1.97 


-0.01 


-1.98 


-2.05 


+0.20 


-1.85 


ZnSe 


-3.30 


-2.85 


-0.38 


-3.23 


-2.05 


-1.17 


-3.22 
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